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O 

' We study baryonic superfluidity under the combined pion condensation with Z\-mixing, 

^ I where both condensates of the neutral and charged pions take place in neutron stars. Outline 

O . of formulation treating such superfluidity is presented and energy gap calculations show the 

' possibility that critical temperatures are higher than about 1 x 10^ K. 

oo . 

CN , (1) Introduction Recent experiments on the spin-isospin excitation mode of 

nuclei^^ have shown that condensations of the neutral pions (tt'^) and the charged 
^ ' pions (vr'^) hkely set in at almost same density not higher than about two times the 

. nuclear density (po = 0.17fm~^) with the assist of the Z\-degrees of freedom. ^-^ This 

view on pion condensation (PC) differs from the previous one that the charged PC 
firstly occurs at about 2po and the neutral PC appears at the higher density, e.g., 
p ~ (3 — 5)po-^'' The combined PC means that both condensates of vr*^ and vr"^ coexist 
in neutron stars (NSs). Such relatively low threshold density of the combined PC 
signifies its existence in most of NSs except for those with very small masses. 

This paper is the first to study the baryonic superfluid (SF) under the combined 
PC with the Z\-degrees of freedom. Here we adopt a model of the combined PC, 
^ ' studied previouly,^^'^^ that both condensates contribute to energy gain without in- 

P5 . terference. The baryon system under this combined PC takes the specific structure 

which we have called the alternating layer spin (ALS) structure. ^^'^^ The SF in such 
a combined PC consists of the pair with the spin-triplet, two-dimensional (2D) P- 
wave and 2D total angular momentum mj = ±2, that is, the pair with the streched 
\ 2D angular momentum. This is because the pair can utilize the attractive effects 

from the LS interaction, needed for the occurrence of the SF at p >P0) as shown 

in previous studies by the present authors. ^^"^^ As this SF contains the ^-P2-pair as 
a dominant one, we call it the ^P2-dominant SF. This paper presents outline of a 
formulation to treat this SF and briefly reports calculated results of energy gaps, 
equivalently critical temperatures. Full description will be reported in a succeeding 
paper. 

(2) Quasi-baryons (n and p) and quasi-particles (rj) in the combined PC The 
ground state of the system without the pairing correlation, denoted as |^^o)i is repre- 
sented by the direct product of the fermion ground state \^f) and the boson ground 
state representing the combined PC: |^o) = I^f) x The \'Pf) is described by 
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the Fock state of the quasi-particles denoted by the ij particles: \<Pf) = 11^°'^'^ 
which tends to the neutron matter without PC. The rj particle signifying the single- 
particle eigen-mode in the combined PC is given by a superposition of a quasi-neutron 
(n) and a quasi-proton (p) constructed in the following way. 

Due to the action of the condensed vr^ field with the momentum ko, the quasi- 
neutron n (the quasi proton p) is formed as a superposition of the neutron n (the 
proton p) and the (A'^) with a same charge and spin component. ^^'^^ Due to 
the action of the condensed field with the momentum k^, a quasi-neutron n (a 
quasi-proton p) is formed as a superposition of the n (p) with the spin= ±1/2 and 
the A- {A++) of the spin state |Z\(±i/2)) = ±k±\A^i/2) /2=f^/3k^\A^y2) /'^■^^ Here 
k± = {kcx i ikcy)/kc = exp{ia9c) for cr = it, where 9c is the angle between kc and 
the X-axis, since wc take the angular-momentum quantization axis along feo in the z 
direction, and kc to be perpendicular to ko, 

The 77-particle is described by a superposition of n and p with opposite spin 
components, because the nucleon-7r"^ interaction flips the spin and isospin. Tranfor- 
mation of the quasi-baryon operators to the r) operators and its inverse are written 
as follows:^) 

Vp = upfip - vppp_, Cp = vphp + uppp_, (la) 

n/3 = uprip + v^Cfi, Pf3- = -Vf^rjis + u^Cp. (lb) 

The (^-states orthogonal to the ry-states are empty because of their energy much 
higher than the ry's. The suffices are given as /3 = {q±,qi^\cr} and /3— = {(qx ~ 
kc),qi^\—cr}, where {q±,qz) denotes the momentum, a/2 is the spin component 
and vr = -|-(— ) means the "^-parity" of the Bloch-orbital state constructed from the 
localized wave functions around layer centers of i =cvcn (odd) numbers. The u, v 
factors satisfy the relation |n^p-|-|u^p = 1. We have the relation, tt = it for a = ±1, 
according to our specification of the ALS state. The coefficients are expicitly given 
as us = u (^) ~ n and vr = —ic^^^^v ~ —ic^^^^v. In the last lines the 

approximation to neglect weak momentum dependence of the u, v factors is used, 
and then u'^ + = 1. In treating the Z\-degrees of freedom in the SF problems, we 
basically follow previous papers.^) 

(3) Quasi-baryon pair operators and pairing interaction Pairing correlation 
under the combined PC works through the side of the Fermi cylinder, namely for 
the {qj_, —q±) pair. Because of the band gap just outside the upper and lower Fermi 
surfaces, the pairing correlation does not work with respect to the momentum q^- 

We consider the quasi-baryon pair operators with the 2D streched angular mo- 
mentum, namely, with the spin ms = a/2 + a/2 = a, the spatial mi = a, and the 
total mj = ms+mL = 2a = ±2 for a = ±1. By this choice we have msniL > which 
gives the attraction from the LS potential because its radial part in the triplet-odd 
state is negative. 

There are two independet groups of the quasi-baryon pairs, the one with mj = 
+2 and the one with mj = —2, because mj is the good quantum number. The 
aspects of the mj = — 2 group are the same with those of the mj = +2 group except 
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that all the components of the 2D angular momenta are reversed. Hereafter we treat 
only the mj = +2 group fixed hy a = +1 for the derivation of the gap equation and 
so forth, although the total energy shift due to the pairing correlation is the sum of 
the contributions from the two groups. 

We define three quasi-baryon pair operators with A = {\ms\ = 1, = 1} and 
mj = +2, for the fin-pair, the pp-pair and the symmetrized (np)-pair, as follows: 

= 4^ / d^g^-^n^ , (2a) 

V2 J V27r Q±,rz ',o- -qj_,-ql ',a 



B^iSi'M^Q.)^^ [ d^,'-^p] (2b) 

V2 J V27r (qx-fec),«i ',o- {-q±-kc),-ql ',a 

+P] (-) (+) I- (2c) 



q±,qz ,0- (-q±-kc),-qz ,<y 



In these B\ J dtpq e^^fi plays a role to make the projection to the 2D orbital angular 
momentum uil = cr, and thus leads to nij = 2a. Owing to this projection operator, 
the quasi-baryon operators are symmetric for — ^ —<lzj and we can restrict to 
> in these pair operators. 
Discriminating these quasi-baryon pair operators by V, we define the most fa- 
vorable pairing interaction in the combined PC as follows: 

q'j_ q± qz V V 

xB'C^;{q'^,q',)B2M,q.). (3) 
By extracting r]^ from B^^^"^ (q±,qz) using Eq. (lb) as 

(+) ^ur/^ (+) , pj (_) ^ -ie^'^'^^v ri^ (_) , (4) 

qj_,qi \<T q^,ql ' ,a {q±-kc),qi ' ,rT qj_,ql ' -a 

we have the r/-pair operators B^^^"^ {q±,qz) expressed by rj^ as follows: 

4rl^(g^, g.) - Bf=;'^ \q^, q,), Bf^(M^,q^ = -v' sfj''^ \q^, q,), and 

B^Z^]{q^,q,)^{-^uv)B^£:f\q^,q,). 

Denoting the coefficients as Cn=+i = u^, Cn=-i = —v^, Cn=o = iuv we have 
the pairing interaction Hamiltonian written in terms of rj and rj^: 

<-a'^ = E E E E E E(^i' r'\c*„'Vi^^^CnW^, r) 

q'j_ 9i q± qz w n 

xbZIH<i'±,i'.)bT2M,1^), (5) 

where the restriction for > is not necessary for the matrix elements of the pairing 
interaction. It is to be noted that the matrix elements are taken with respect to the 
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states of the quasi-baryon pair, not to those of the ij pair, despite of the one to one 
correspondence between V and 11. 

(4) Bogoliubov transformation and coupled gap equation In treating a coupled 
pairing of three channels with respect to 11, we apply the Bogoliubov transformation 
following the method previously developed in the study of the ^P2-pairing, referred 
to as T70.^°) Namely, the superfluid groud state \^o) is obtained by operating a 
unitary tansformation e*'^ on the nonsuperfluid ground state \^f)- 

|<Fo) = e'^l^p) = e'^^"'J=+^ X e^'^^'"J=-2|#F). (6) 

We first express iSAmj=+2 part using B^^^^ and B^^^ and then rewrite it in 
terms of r) and r/^, using q = {qj_, q^), as follows: 

(7) 

The suffix a can be suppressed because cr = +1 is fixed for mj = +2, and the 
transformation is carried out with respect to vr. As for the formal proccdTirc, we 
can proceed in pararell with T70, treating vr in place of a. The 9 function has the 
symmetric and parity-odd property: 6{q, 7^1,1^2) = 9{q, 1:2,^1) and 6'(— g, tti, 7r2) = 
— 0(qf, TTi, 7r2). We write the operators rj and rj^ columnwise as 

and similarly write the transformed operators a and , which we abbreviate as 
SFQP for the superfluid quasi-particle. Then, denoting ^(g,7ri,7r2) by the (2 x 2) 
matrix 0{q), we have the Bogoliubov transformation: 

dq = e'^^-^^qC-'^^-^ = U{q)ffq - V{q)fi^q ^ (9a) 

dq t = c'^-^r^j ^Q-iSAmj = U{q)ffq t _ y t (q)fj_g. (9b) 
U and V are the (2 x 2) matrices with respect to vr, which are given by 

U{q) = cos6)'(q), (10a) 



V{q)^0'{q)-'sm0'{q)0{q), (10b) 

where 0'{q) is the Hermite matrix determined by the relation 6>'(g)^ = 0{q)0{q) 
We rewrite the pairing interaction Hamiltonian using the inverse Bogoliubov 
transformation in terms of a and and lead it into the reduced form consisting of the 
constant term without and a, the single SFQP term with a^Q, and the dangerous 
term involving aa and a^a^ . From the vanishing of the dangerous term, we reach a 
coupled gap equation among three gap functions An{q±.,qz) of il = +1, —1, 0: 
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where E{q±,qz) = ^(^,{qL.,qz)+D'^{q^,qz), l,^{q±,qz) ^ fi^q'j_/2Mnm* (m* being 

an effective mass parameter of the ?7-particle), Kr[=^i^-i^o = 1, 1,2 and the squared 
energy gap D'^ = Kw^'n/'iTr. In getting this familiar form of the SFQP energy, 
we have taken an approximation that a squared gap matrix {D"^) is replaced by an 
appropriate c-number form for which we choose its trace. 

(5) Estimate of energy gap and critical temparature For the NS matter in the 
combined PC, the n is the dominant component, even though the p is generated by 
the charged PC, and in the pairing problem the IJ = +l-channel corresponding to 
the nn-pair contributes most importantly. Generally, small components play a role 
to enhance the main component in the coupled gap equation. Typical examples can 
be seen in the tensor-force coupling in the + ^-F2-pairing^^) and the ^Si + ^Di- 
pairing.^^ Here the 77 = —1 and 77 = channels enhance the gap function of the 
dominant 77 = +1 channel through the coupling terms. Thus we can estimate the 
lower limit of energy gaps by making calculations in the 77 = +1 channel only. 

The matrix elements in the main channel {V = nh) is approximately given by 
the following 2D radial integral, after properly managing the normalization factors, 
as follows: With use of the cylindrical coordinate {r±,(p,z), 

{Q±,Q.,r\C!V^i"''^\q'^,q:,V) ^^^^ H dr^r^J,{q^r^)4''"''\r^)MqW), 

+ Vi) Jo 

(12) 

where v'^ is the p-mixing probability, yf the A -mixing probability and Ji{q_Lrj_) 
the Bessel function. Here we take an approximation suitable for the well-developed 
ALS structure, because this structure grows just beyond the on-set density of the 
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Fig. 1. The thick sohd curve shows the efltec- 
tive 2D pairing potential of the main hn 
channel in the combined PC obtained for 
OPEG-B potential and at p = 3.5po- The 
corresponding potential in the neutral PC 
(the one multiplied by the attenuation fac- 
tor=0.500) is shown by the solid thin (dot- 
ted) curve. The dash-dotted curve shows 
the attraction due to the Zi~-mixing 
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Fig. 2. The solid (dotted) curves show the 
density dependence of SF critical temper- 
atures Tc in the combined PC (the neu- 
tral PC) for OPEG-B potential. The num- 
bers attached to the curves are the effective 
mass parameters mJJ. 
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neutral PC.^^ Then the qz dependence becomes irrelevant. (rj,) is the 2D 

effective potential obtained after the integation over z and 0, where we suppress the 
q± dependence although include it in solving the gap equation. In the single channel 
we have D{q±) = A-^-l{q±) / ^/Ttt . 

Energy gap calculations are done by adopting the baryon-baryon potential work- 
ing in the {N + A) space extended from NN OPEG-B potential^^^ and by using the 
PC parameters for the strong PC case given in refs. 3), 13) but slightly modified so 
as to have smooth density dependence. Fig.l shows an example of the 2D effective 
potential calculated at p = 3.5po- Attenuation of the pairing attraction caused due 
to the charged PC is severe as (1 - + yf)"^ ^ 0.74 0.50 0.63, for such 

density variation as ~ 0.10 0.24 0.17 and yj ~ 0.045 ^ 0.072 0.042 at 
/9 = (2 ^ 4 ^ 6)po- But the Z\~-mixing brings about the short-range attraction 
shown by the dot-dashed curve in Fig. 1, and the pairing attraction is recovered 
from the attenuated one considerably as illustrated by the thick arrow. As a result, 
the reduction of the pairing attraction from that in the neutral PC is gentle. 

Superfluid critical temperatures Tc obtained from energy gap calculations, 
ksTc — 0.57D(q±F) {q±F being the Fermi momentum), are shown in Fig. 2 at 
densities p = {2 — 6)po and for several values of the effective mass parameters of 
the r?-particle m*. Since the m* in the combined PC are about (0.8-0.9), we have 
Tc >1 X 10^ K for this interaction. The result is meanigful, because such Tc efficiently 
suppress too rapid cooling due to the pion direct Urea process taking place in most 
of NSs, except for those with very small masses. This point is to be considered in 
the studies of NS cooling. ^^^^^^^ To proceed further, wc need detailed studies on the 
matrix elements of the channel coupling and solutions of the coupled gap equation, 
which will be reported elsewhere. 

The authors thank S. Tsuruta and T. Tatsumi for their cooperative discussions. 
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